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PREFACE 


The  research  reported  in  the  paper  reprinted  within  was  started  by  the  two 
authors  at  the  U.S.  Naval  Research  Laboratory,  Washington,  D.C.  Soon 
after,  one  of  the  authors  (W.A.  Kuperman)  transferred  to  SACLANTCEN  where 
he  performed  most  of  his  work.  In  addition,  a  complete  computer  model  for 
performing  the  calculations  has  been  developed  at  SACLANTCEN. 
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Spatial  correlation  of  surface- generated  noise  in  a  stratified 
ocean 

W.  A.  Kuperman 
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and  Naval  Research  Laboratory,  Washington,  D.C  20375 


F.  Ingenito 

Naval  Research  lAtttftpt 
(Receiv«T$May  1979; ; 


try,  Washington,  D.C.  20375 
Iccepted  for  publication  3  March  1980) 


A  model  is  developed  for  the  calculation  of  the  spatial  properties  of  the  noise  field  produced  in  a 
stratified  ocean  by  the  action  of  wind  at  the  surface.  The  random  noise  sources  are  represented  by 
correlated  monopoles  distributed  over  an  infinite  plane  located  an  arbitrary  depth  below  the  surface. 
Wave-theoretical  methods  are  applied  to  derive  expressions  for  the  intensity  and  spatial  correlation  of  the 
noise  field.  A  normal-mode  representation  of  the  noise  field  is  used  to  reduce  these  expressions  to  forms 
which  allow  physical  interpretation  and  are  suitable  for  numerical  computation.  Examples  are  given  of 
intensity  profiles  and  spatial  correlation  in  the  vertical  for  three  generic  sound-speed  profiles.  The  results 
show  that  the  sound-speed  profile  and  the  presence  of  the  bottom  can  be  important  in  determining  the 
spatial  properties  of  the  noise  field.  An  example  is  given  of  a  calculation  of  the  horizontal  spatial 
correlation  using  the  fast  field  program  (FFP). 

PACS  numbers:  43.30.Nb,  43.30.Cq 


INTRODUCTION 

Detection  of  acoustic  signals  in  the  ocean  is  always 
performed  against  a  noise  background.  Arrays  of  sen¬ 
sors  may  provide  some  discrimination  against  noise, 
the  degree  of  discrimination  being  expressed  by  the  ar¬ 
ray  gain.  The  array  gain,  defined  as  the  ratio  of  the 
signal  to  noise  of  the  array  output  to  the  signal  to  noise 
of  the  output  of  a  single  element,  can  be  shown  to  de¬ 
pend  on  the  spatial  correlation  of  the  noise  field.1  in  ad¬ 
dition,  recently  developed  optimal  array  processing 
techniques2*3  require  knowledge  of  the  spatial  correla¬ 
tion  of  the  noise  field. 

One  of  the  major  components  of  the  ambient  noise 
field  in  the  ocean  is  produced  by  the  action  of  the  wind 
at  the  surface.  Previous  theoretical  studies  of  the  spat¬ 
ial  structure  of  surface  generated  noise  have  been  car¬ 
ried  out  with  deep  water  applications  in  mind.4*5  Thus, 
the  ocean  has  been  modeled  as  a  homogeneous  half¬ 
space  which  allows  straight-line  propagation  without  re¬ 
flection,  greatly  simplifying  the  calculation.  Such  mod¬ 
els  are  of  doubtful  validity  in  shallow  water  where  the 
acoustic  field  interacts  strongly  with  the  bottom.  Cox® 
has  pointed  out  that  the  assumption  that  noise  arrives 
only  from  above  the  horizontal  is  counter  to  experi¬ 
mental  evidence.  He  shows  how  the  spatial  structure  of 
the  noise  field  is  related  to  a  plane-wave  directivity 
function  in  terms  of  a  sum  of  angular  harmonics.  The 
coefficients  of  the  harmonics  can  in  turn  be  related  to 
deep  water  experimental  results.  We  note  however  that 
in  deep  or  shallow  water  the  sound  speed  is  not  constant 
in  depth,  a  fact  which  may  have  a  profound  effect  on  the 
noise  field,  as  it  does  on  the  signal  field.  Since  the 
ocean  is  horizontally  stratified  it  is  quite  possible  that 
the  acoustic  field  cannot  be  expressed  in  terms  of  the 
same  weighted  set  of  plane  waves  (same  directivity 
function)  at  each  of  the  hydrophones  of  an  array,  partic- 
liarly  a  large  aperture  vertical  array. 
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In  this  paper,  using  wave  theory,  we  develop  a  model 
of  surface  generated  noise  in  which  the  ocean  is  strati¬ 
fied  in  depth.  The  acoustic  properties  of  the  ocean  bot¬ 
tom  are  included,  as  are  the  statistical  properties  of 
the  surface.  In  Sec.  I  we  derive  a  general  expression 
based  on  wave  theory  for  the  cross-spectral  density  of 
surface  generated  noise.  In  Sec.  II  we  apply  this  form¬ 
alism  using  a  normal- mode  representation  of  the 
Green’s  function  for  the  problem.  This  allows  us  to 
gain  some  physical  insight  into  how  the  noise  is  spatial¬ 
ly  distributed.  Section  HI  presents  some  numerical  re¬ 
sults  for  realistic  ocean  environments.  In  Appendix  A 
we  show  that  the  results  derived  in  Sec.  I  reduce  to 
earlier  results4  when  the  appropriate  limits  are  taken. 
Finally,  Appendix  B  presents  some  purely  analytic  re¬ 
sults  for  an  idealize’  waveguide. 

I.  DERIVATION  OF  THE  SPATIAL  PROPERTIES  OF 
SURFACE  GENERATED  NOISE 

The  model  geometry  is  shown  in  Fig.  1.  The  figure 
illustrates  a  simple  case,  with  a  layer  of  water  overly¬ 
ing  a  semi-infinite  bottom,  the  density  and  sound  speed 
of  the  water  and  the  bottom  given  by  p,,  c, (z)  and 
pj.cjte),  respectively.  The  theory  to  be  presented  be¬ 
low  is  also  applicable  to  more  complex  environments. 
For  example,  attenuation  in  the  water  and  the  bottom 
can  be  included  as  can  a  layered  bottom  with  finite  rig¬ 
idity.  Figure  1  is  merely  meant  to  suggest  that  the  en¬ 
vironment  must  be  stratified  in  depth,  thus  ensuring 
separability  of  the  wave  equation.  In  the  development 
below,  we  drop  the  subscripts  distinguishing  the  sound 
speeds  in  the  water  and  bottom  and  denote  the  sound 
speed  anywhere  in  the  medium  oy  c(z). 

Consider  an  infinite  plane  parallel  to  the  surface  and 
located  below  the  surface  at  depth  z'.  Assume  that  at 
each  point  in  the  plane  there  is  a  monopole  source  of 
strength  six',  /),  where  r'  is  the  radial  vector  in  the 
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FIG.  1.  The  model  geometry  showing  the  source  plane,  at 
depth  z‘  below  the  surface,  and  the  two  field  points  (rj.Zj)  and 
(fj.Zjh 


source  plane  and  t  is  the  time  variable.  Let  the  func¬ 
tion  s(r',t)  be  a  random  variable.  These  monopoles  will 
couple  into  the  water  column  as  dipoles  because  of  the 
pressure  release  surface;  this  effect  is  automatically 
incorporated  in  the  wave-theoretic  treatment  used.  We 
use  monopole  sources  because  they  represent  the  basic 
fluctuating  volume  source7  and  more  complicated 
sources  can  be  considered  to  be  a  sum  of  these  sources 
appropriately  distributed  in  space.  Therefore,  the 
source  function  is  s(r',f)C(z-z')  so  that  the  field  (vel¬ 
ocity  potential)  In  the  water  column  *(r,z,  f)  satisfies 


the  wave  equation 

where  ate)  is  the  Dirac  delta  function. 

We  represent  *  and  s  by  their  Fourier  transforms 

#(r,z,f)  =  (2i)*1/*  f  doJip,Jr,z)exp(-iwt) ,  (2) 

s(r',f)  =  (2ir)'l/*  jT  da»su(r')exp(-ifc>ri,  (3) 

where  u>  denotes  angular  frequency. 

Inserting  Eqs,  (2)  and  (3)  intoBq.  (1),  we  obtain,  after 
some  manipulation, 

(v* +*»)?„ =-Sw(r')6(z-z');  fe=«/c(z).  (4) 

Equation  (4)  has  the  solution 

<pjr,z)~  f  dxT'Sjr')G{T,T,;z,z‘),  (5) 


where  G(r,r';z, «'),  the  Green’s  function  of  the  prob¬ 
lem,  satisfies  the  Helmholtz  equation 

(V*  +**)G(r,  r';  z,  z')= -(l/r)6*(r  -  r')6(z  -  « ')  (6) 

and  the  appropriate  boundary  conditions.  Equation  (S) 
simply  states  that  the  total  velocity  potential  is  obtained 
by  summing  over  all  source  contributions.  We  note 
here  that  Su  is  the  spectral  strength  of  the  noise 
sources  and  that  the  total  field  is  given  by  integrating 
over  all  frequencies  as  stated  in  Eq.  (2).  In  order  to 


simplify  notation,  we  shall  usually  drop  the  subscript 

< i j. 

The  cross-spectral  density  is  a  measure  of  the  spatial 
coherence  of  the  noise  field.  To  obtain  the  cross-spec¬ 
tral  density  we  form  the  product  of  and 

</>*{ra,  za)  and  take  the  ensemble  average  (tp*  is  the  com¬ 
plex  conjugate  of  <p).  Thus, 

( ^(r,,z,V*(ra,  za)) 

=  /  /  ^V*r»<S(r')S*(r")) 

xG(r1,r';z1,z')G*(ra,r";za,z'),  (7) 

where  the  angle  brackets  indicate  an  average  taken  over 
the  random  function  S.  It  will  be  convenient  to  use  a 
transverse  Fourier  representation  of  the  Green’s  func¬ 
tion®  which  we  write  as 

G(r,r';z,z') 

=  -2~-  J  d2»jg(ii;a,a')exp(i»)-(r-r')] ,  (8) 

where  g(jj;z,z')  satisfies  the  equation 

+  \k2(z)-ri2]g=--^  6{z-z'),  (9) 

which  follows  from  Eq.  (6). 

Using  this  Fourier  representation,  we  can  express  the 
cross-spectral  density  function  of  the  noise  field  as 

( (t j ,  Zj)^>®(ra,  Z2)) 

=  f f  d*r' d*r0(S(r')S(r*)) 

X  1(2*)*  1 1 d **  d*n  '8^r,;  lx,t  )g*  ^  ’ **’ z 

*  expjzij  •  (rj  -  r')]  exp[-zij'  •  (r2  -  r')] .  (10) 

Now  let  R  =r,  - ra  and  p=r'-  r",  and  assume  that  the 
spatial  coherence  of  the  noise  sources,  (S(r')S(r")>,  de¬ 
pends  only  on  p.  We  denote  <S(r')S(r")>  as  qxN(p).  Sub¬ 
stituting  for  r,  and  r'  in  Eq.  (10),  the  integrations  over 
r"  and  ij'  can  be  performed,  resulting  in 

CW(R, Z},  z2)  *  (*p(rj,  Zi)*f>*(T2t  z2)) 

~q2 1 J  d2PdxtiN(p)g[i);zltz') 

x£*(»?;*2»*')exp|*s  *  (R  —  p)] ,  (ii) 

where  w  is  used  as  a  subscript  to  remind  us  that  the 
cross-spectral  density  function  depends  on  frequency. 
Since  g  and  g*  depend  on  the  magnitude  of  jj,  but  not  its 
direction  (see  Eq.  (9)],  we  can  perform  the  integration 
over  the  azimuthal  angle  associated  with  tj,  with  the  re¬ 
sult  that  the  cross-spectral  density  function  takes  the 
form 

CU(R,  zl,z2)  =  2zqx  J dpN(p)  J  T) 

X£*(n:z.,c')J0<T?|R-p!),  (12) 

where  J0  is  the  Bessel  function  of  zero  order.  Another 
form  for  C„,  which  is  particularly  simple,  can  be  ob¬ 
tained  by  expressing  N'  p)  by  its  Fourier  transform 
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Pbl): 

N(p)=£  f  dhiPWe1''*  •  U3) 

Then  the  integration  over  p  in  Eq.  (11)  can  be  per¬ 
formed,  giving 

CU(R,  zlt  z2)=2vti2  j d2ri P{i))g(t];  z ') 

*  g*(V-,z2.z‘ ')«'»'*  .  (14) 

Equation  (12)  can  be  put  into  a  form  which  will  be  useful 
for  later  calculations  by  decomposing  the  Bessel  func¬ 
tion  into  a  sum  of  Hankel  functions 

J0U)  =  M1>(z)+H'2>(z)),  (15) 

where  the  superscripts  denote  the  Hankel  function  of 
first  and  second  kind.  Using  the  relation  -H^K-z) 
={i£\z)  a*:d  noting  that  g  and  g*  are  even  in  i)  we  can  ex¬ 
tend  the  integration  over  ;;  from  -»  to  ».  Equation  (12) 
becomes 

CW(R,  zt,z2)  =  vq2  j  dp  Nip)  J  rjdtjH^KnlR-p\) 

xglT/iz^z'^tozj.z').  (16) 

From  the  expressions  for  the  cross-spectral  density 
function  of  the  noise  field  [Eqs.  (12),  (14),  and  (16)], 
we  can  immediately  make  some  comments  about  its 
structure.  In  any  horizontal  plane,  it  is  independent  of 
absolute  position  and  depends  only  on  the  horizontal 
vector  R  connecting  the  field  points.  In  the  vertical,  the 
spatial  coherence  depends  not  only  on  separation  dis¬ 
tance  but  also  on  the  absolute  depth  of  the  field  points. 
Hence,  in  general,  the  noise  is  not  spatially  stationary 
in  the  vertical. 

An  important  special  case  is  that  of  uncorrelated 
noise  sources.  Equation  (16)  is  an  expression  for  the 
cross- spectral  density  function  of  the  noise  field  as  a 
function  of  the  spatial  coherence  of  the  noise  sources 
N(p).  For  uncorrelated  noise  sources,  it  has  been 
shown  that® 

Mp)  =  26{p)/k2p.  (17) 

Using  Eq.  (17)  in  Eq.  (16)  we  get 

CuQR,zl,zz)=4z2q2k~2  f  q  dqH^’iqR) 


enclature  that  the  correlation  function  of  the  noise  field 
is  given  by 

(4>(r,,z„  l)4-*(r2,22,  /  +  t)) 

=  J  Q,(R,zl,z2)exp[-iu>T)dw ,  (21) 

«e 

where  t  is  a  time  delay. 

In  this  paper,  we  are  mainly  concerned  with  the  dis¬ 
tribution  of  noise  in  a  stratified  ocean,  and  in  particular 
in  situations  where  the  acoustic  properties  of  the  ocean 
bottom  have  a  profound  effect  on  the  acoustic  field. 
Nevertheless,  the  model  should  also  handle  situations 
where  the  bottom  is  not  important,  for  example,  the 
deep  ocean.  In  Appendix  A  we  show  that  the  above  theo¬ 
retical  results  reduce  to  earlier  work®  where  the  ocean 
was  modeled  as  a  semi- infinite  isovelocity  half-space. 

II.  NORMAL  MODE  REPRESENTATION  OF  THE 
NOISE  FIELD 

In  this  section  we  apply  the  results  of  Sec.  I  to  a 
stratified  medium,  that  is  a  medium  in  which  the  sound 
velocity  and  dei.sity  of  the  medium  are  functions  of 
depth  z  only.  The  Green’s  function  can  be  expressed  in 
several  equivalent  ways;  in  this  section  we  use  a  nor¬ 
mal-mode  representation  in  which  the  Green’s  function 
is  expanded  in  terms  of  the  normal  modes  of  the  sys¬ 
tem.  If  the  medium  is  finite  in  depth  with  appropriate 
conditions  given  at  the  boundary,  the  normal  modes  will 
be  discrete  and  the  propagating  modes  will  be  finite  in 
number.  However,  if  the  medium  is  infinite  in  depth, 
there  will,  in  general,  exist  a  finite  number  of  discrete 
modes  and  an  infinite  set  of  continuous  modes.  The 
Green’s  function  expansion  will  then  consist  of  a  dis¬ 
crete  sum  plus  an  integral  over  the  continuous  modes. 

For  simplicity,  we  restrict  ourselves  to  that  part  of 
the  noise  field  which  can  be  represented  by  a  discrete 
set  of  normal  modes.  From  the  above  discussion  we 
see  that  this  will  be  a  complete  description  for  the 
pressure- release/ rigid  waveguide  discussed  in  Appen¬ 
dix  B,  but  not  for  more  realistic  ocean  models.  The 
latter  usually  consist  of  a  layer  of  water  and  several 
sedimentary  layers  overlying  a  semi-infinite  basement. 
However,  by  making  the  acoustic  impedance  of  the 
basement  very  high,  we  can  minimize  the  importance 
of  the  continuous  modes. 


x  git);z„z')g*in;z2,z').  (18) 

When  we  set  R  =0  and  zi=z1=e  in  the  expression  for 
Cu,  we  obtain  a  quantity  proportional  to  the  intensity  of 
the  noise  field  at  a  point.  Equation  (16)  then  reduces  to 


The  Green’s  function gii)‘,z, z')  can  be  written  in  terms 
of  the  normal  modes  as  follows®: 


giT]-,z,z')  = 


pJz  2  y 

2»  V 


UJz')Um(z) 


(22) 


C,j{0,z,z)siu{z)  =  vq3  j dpN(p) 

X  ^"r)dr]Hi,,(ijp)lg(jj;zfz')l2  (19) 

or  alternatively,  from  Eq.  (14), 

j d2ii  P(^)lgiij;z,z')l2 .  (20) 

The  expressions  given  by  Eqs.  (18),  (19),  and  (20)  will 
be  useful  later.  We  also  mention  for  clarity  in  nom- 


where  VJ(z)  and  k„  are  the  normalized  mode  amplitude 
function  and  the  wavenumber  of  the  «th  mode  and  are 
solutions  of  the  eigenvalue  problem  defined  by  the  equa¬ 
tion 

^^&  +  \k2(z)-k2Mz)=0,  (23) 

with  the  appropriate  boundary  conditions.  In  Eq.  (22) 
ps(z ')  is  the  density  of  the  medium  at  the  depth  z'  and 
k(z)  =  u/ciz)  with  clz)  being  the  sound  speed. 
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We  assume  that  k„  is  a  complex  number  of  the  form 

km  =  *,,  +  iam,  (24) 

with  k„,  a„  fe  0;  an,  the  imaginary  part  of  k„,  is  the 
modal  attenuation  coefficient.  It  is  interesting  to  note 
that  we  must  include  attenuation  in  the  system  to  obtain 
a  finite  cross-spectral  density  function.  This  is  be¬ 
cause  sound  trapped  by  the  layered  medium  (represent¬ 
ed  by  the  discrete  modes)  suffers  cylindrical  spreading 
while  the  amount  of  energy  radiated  by  the  noise 
sources  increases  as  the  square  of  the  range  from  the 
field  points.  Hence,  the  contribution  to  the  intensity  of 
distant  sources  increases  with  range  and  the  total  in¬ 
tensity  diverges.  Any  amount  of  attenuation  in  the  sys¬ 
tem  will  cause  the  intensity  to  decay  exponentially  with 
range  and  ensure  convergence.  It  is  important  to  note 
that  the  resulting  cross-spectral  density  functions  and 
intensities  will  depend  on  the  attenuation  chosen. 

We  now  insert  the  Green’s  function  of  Eq.  (22)  into  the 
expression  for  the  cross-spectral  density  function 
[Eq.  (16)]  and  evaluate  the  tj  integral.  From  Eq.  (22) 
we  see  that  the  integral  of  Eq.  (16)  has  simple  poles  at 

tj  =  ±k„  ±k* ,  (25) 

of  which  the  poles  at  and  -Jt£  are  in  the  upper  half¬ 
plane.  Using  standard  methods  of  complex  integration 
we  close  the  contour  in  the  upper  half-plane  with  a 
semicircle  of  large  radius  and  evaluate  the  residues  of 
the  integral.  The  result  is 

Cuq».  z,)  =■&&<*')  fd(tN(p) 

*T,V«U’)Um(zl)UJz‘)UmU,) 


(26; 

where 

/•.  =  !/(** -***)•  (27) 

The  quantity  /_  is  a  measure  of  the  coherence  be¬ 
tween  the  normal  modes  which  make  up  the  noise  field. 
For  example,  if  /„,  vanishes  for  n  *  m,  then  the  noise 
field  reduces  to  an  incoherent  sum  over  the  normal 
modes.  Writing  /_  in  terms  of  the  complex  fe.’s  and 
assuming  that  *„»  o„,  *m  »am,  we  get 

{1/U2.-k.)  for  m*n, 

,/a-  ,  (28) 

l/4tayt,  for  h  =  m . 

We  see  that  tire  n  =  m  terms  in  the  sum  of  Eq.  (26)  be¬ 
come  infinite  in  the  absence  of  attenuation.  This  is  due 
to  the  contributions  of  distant  sources,  as  discussed 
above.  The  n*m  terms  remain  finite  because  they  are 
products  of  different  modes  with  rapidly  oscillating 
phases  which  give  negligible  contribution  from  distant 
sources.  Equation  (28)  also  indicates  that  if  the  attenua¬ 
tion  coefficients  a„  are  much  smaller  than  the  smallest 
separation  between  eigenvalues,  then  the  noise  field  can 
be  approximated  by  an  incoherent  sum  of  modes.  This 
is  often  the  case  in  shallow  water.  We  can  further  sim¬ 
plify  Eq.  (26)  by  approximating  km  by  its  real  part 


Then  Eq.  (26)  reduces  to 


Cy(R| 


IpMI 


j  dpN(p) 


(29) 


where  we  have  neglected  the  n*m  terms.  Finally,  when 
the  noise  sources  are  completely  uncorrelated,  the 
cross-spectral  density  takes  the  simple  form 

=  E  W  J^ft) .  (30) 

From  Eq.  (29)  and  Eq.  (30)  it  is  obvious  that  the 
structure  of  the  noise  field  is  highly  dependent  on  the 
attenuation;  in  shallow  water  the  attenuation  is  usually 
dominated  by  the  acoustic  interaction  with  the  bottom 
sediments. 

In  Appendix  B  we  evaluate  the  cross-spectral  density 
function  for  a  case  which  can  be  done  analytically:  an 
isovelocity  waveguide  bounded  above  by  a  pressure- re¬ 
lease  surface  and  below  by  a  rigid  bottom.  Though  it  is 
not  very  descriptive  of  a  real  ocean  environment,  the 
analytic  calculations  are  helpful  in  understanding  how 
surface  noise  is  distributed  in  a  waveguide. 


III.  NUMERICAL  RESULTS  AND  EXAMPLES 

In  this  section  we  present  sample  calculations  which 
exhibit  some  of  the  properties  of  the  spatial  correlation 
and  intensity  of  the  noise  field.  Most  of  the  calculations 
were  made  using  the  normal-mode  representation  of  the 
noise  field  presented  in  Sec.  II,  but  we  emphasize  that 
the  model  presented  here  is  not  bound  to  a  specific  rep¬ 
resentation;  any  wave-theoretical  representation  can  be 
used.  As  an  example,  we  will  also  give  some  results 
calculated  by  a  modification  of  the  fast  field  program10 
(FFP),  where  the  function  g(ij\z,z')  is  calculated  di¬ 
rectly  and  Eqs.  (18)  and  (19)  are  used. 

First  we  present  three  cases  which  illustrate  the  ef¬ 
fect  of  sound-speed  profile  and  frequency  on  the  intens¬ 
ity  and  spatial  coherence  of  the  noise  field.  In  all  three 
cases  the  water  depth  is  50  m  and  the  bottom  consists 
of  20  m  of  sediment  overlying  a  hard  basement.  The 
sound  speed,  density,  and  attenuation  coefficient  of  the 
sedimentary  layer  are  characteristic  of  sand-silt- 
clay."  The  noise  sources  are  assumed  to  be  uncorre¬ 
lated  and  located  0.5  m  below  the  surface;  they  are 
equivalent  to  sources  at  the  surface  with  cose  direction¬ 
ality. 

The  three  sound-speed  profiles:  isovelocity,  down¬ 
ward  refracting,  and  upward  refracting  are  shown  in 
Fig.  2.  The  corresponding  noise  intensities  are  shown 
in  Figs.  3,  4,  and  5  as  functions  of  depth  for  the  fre¬ 
quencies  200,  400,  and  800  Hz.  The  results  for  the  iso¬ 
velocity  and  downward-refracting  cases  are  qualiatively 
similar.  Both  cases  show  decreasing  noise  intensity  as 
a  function  of  depth,  with  a  faster  rate  of  decrease  for 
the  higher  frequencies.  The  intensity  decrease  is 
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FIG.  2.  Isovelocity,  downward-refracting,  and  upward-re¬ 
fracting  sound-speed  profiles  used  in  the  calculations. 


caused  by  the  frequency-dependent  attenuation  of  the 
system.  For  the  down  ward- refracting  profile  the  atten¬ 
uation,  which  is  dominated  by  the  bottom,  is  greater, 
resulting  in  a  more  rapid  decrease.  The  peak  in  intens¬ 
ity  which  appears  at  all  three  frequencies  in  the  up¬ 
ward-refracting  case  is  caused  by  the  low-order  modes, 
which  are  the  dominant  contributors  to  the  intensity. 

The  low-order  modes  are  trapped  in  the  upper  part  of 
the  water  column  and  hardly  interact  with  the  bottom. 
Thus,  their  attenuations  are  very  &mall  and  contribu¬ 
tions  from  distant  sources  are  important.  The  domin¬ 
ant  low-order  modes  are  strongest  in  the  upper  part  of 
the  water  column,  resulting  in  the  observed  peak. 

For  the  noise  intensity  plots  shown  in  Figs.  3,  4,  and 
5  no  absolute  levels  are  given.  The  model  does  not  pre¬ 
dict  the  levels  of  the  noise  sources,  expressed  by  q2  in 
Eq.  (11),  which  we  expect  to  be  dependent  on  frequency. 
For  the  purposes  of  these  calculations  q 2  has  been  set 
equal  to  unity  in  all  cases. 

Figures  6,  7,  and  8  show  the  spatial  correlation  func¬ 
tion  for  the  same  cases  as  above.  (For  a  single  fre¬ 
quency  and  zero  time  delay  the  spatial  correlation  func¬ 
tion  is  equal  to  the  real  part  of  the  cross-spectral  dens¬ 
ity.)  For  comparison  the  results  for  surface  sources 
having  cosp  directionality  in  a  semi-infinite  homogen¬ 
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FIG.  4.  Noise  intensity  as  a  function  of  depth  for  the  down¬ 
ward-refracting  profile  shown  in  Fig.  2(b)  and  for  the  frequen¬ 
cies  200,  400,  and  800  Hz. 


eous  medium  (Cron  and  Sherman4  and  Appendix  A)  are 
also  shown.  Again  '  le  upward- refracting  case  is  the 
most  interesting,  i.i  this  case,  as  mentioned  above,  a 
few  low-order  modes  dominate  the  noise  field,  result¬ 
ing  in  high  values  of  coherence  throughout  the  water 
column.  In  the  isovelocity  and  downward- refracting 
cases  many  modes  contribute  to  the  noise  field  and  the 
coherence  is  much  closer  to  the  Cron  and  Sherman  re¬ 
sults. 

Finally,  we  give  a  calculation  of  the  spatial  correla¬ 
tion  in  the  horizontal  direction.  We  have  assumed  an 
isovelocity  water  layer  100  m  thick  with  a  sound  speed 
of  1500  m/s  and  a  single  semi-infinite  bottom  of  unit 
density,  sound  speed  of  1600  m/s,  and  an  attenuation 
coefficient  of  1  dB/A.  Figure  9  shows  the  horizontal 
spatial  correlation  at  100  Hz  along  with  the  Cron  and 
Sherman  result  for  comparison.  For  this  environment 
the  model  spectrum  consists  of  a  discrete  part  and  a 
continuous  part,  both  of  which  contribute  to  the  noise 
field.  The  correlation  was  calculated  using  a  combina¬ 
tion  of  the  normal-model  and  FFP  methods.  Thus.  Eq. 
(26)  was  used  for  the  discrete  normal-mode  part  and  a 
modification  of  Eq.  (18)  was  used  for  the  continuous¬ 
mode  part.  In  Fig.  10  we  have  plotted  the  discrete  and 
the  continuous  contributions  separately,  both  normal¬ 
ized,  to  illustrate  the  differences  between  the  two  con¬ 
tributions.  The  continuous  part,  while  more  coherent 
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FIG.  6.  Vertical  spatial  correlation  functions  for  the  isovelo¬ 
city  profile  shown  in  Fig.  2(a)  as  a  function  of  DA,  where  D  is 
foe  receiver  separation,  X  the  acoustic  wavelength,  and  with 
one  receiver  fixed  at  20-m  depth.  Three  frequencies  are 
shown:  —  200  Hz;  •••  400  Hz;  and  800  Hz.  Also 

shown  is  foe  result  for  a  semi-infinite  homogeneous  medium 
calculated  from  Eq.  (A22): - 


spatial  correlation 

FIG.  7.  Vertical  spatial  correlation  functions  for  foe  down¬ 
ward-refracting  profile  shown  in  Fig.  2(b)  as  a  function  of  DA. 
where  D  is  foe  receiver  separation,  X  the  acoustic  wavelength, 
and  with  one  receiver  fixed  at  40  m.  Three  frequencies  are 
shown:  —  200  Hz;  •••  400  Hz;  and  800  Hz.  Also 

shown  Is  foe  resuU  for  a  semi-infinite  homogeneous  medium 
calculated  from  Eq.  (A22):  — . 
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FIG.  8.  Vertical  spatial  correlation  functions  for  foe  upward- 
refracting  profile  shown  in  Fig.  2(c)  as  a  function  of  DA, 
where  D  is  foe  receiver  separation,  X  foe  acoustic  wavelength, 
and  with  one  receiver  fixed  at  10  m.  Three  frequencies  are 
shown:  —  200  Hz;  •  - •  400  Hz;  and  800  Hz.  Also 

shown  is  foe  result  for  a  semi-infinite  homogeneous  medium 
calculated  from  Eq.  (A22): - . 

for  small  receiver  separations  (relative  to  a  wave¬ 
length),  quickly  becomes  less  coherent  than  the  discrete 
part,  the  latter  maintaining  some  degree  of  coherence 
over  several  wavelengths. 

The  relative  importance  of  the  discrete  and  continuous 
parts  of  the  normal-mode  spectrum  is  dependent  on  the 
total  loss  of  the  system.  For  low  loss  the  discrete 


FIG.  9.  Horizontal  correlation  function  (solid  line)  for  an  iso¬ 
velocity  water  layer  overlying  an  isovelocity  semi-infinite  bot¬ 
tom  as  a  function  of  DA,  where  D  is  foe  receiver  separation 
and  X  is  foe  acoustic  wavelength.  The  dashed  tine  is  foe  result 
for  a  semi-infinite,  homogeneous  medium  (Eq.  (A18)  with  m 
=  ll. 
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FIG.  10.  The  horizontal  spatial  correlation  for  the  same  case 
as  Fig.  9  showing  the  discrete  ( — )  and  continuous  ( - )  con¬ 

tributions.  both  normalized. 

modes  dominate,  as  they  can  be  propagated  very  large 
distances  from  a  very  large  area.  In  high- loss  cases 
the  continuous  modes  tend  to  dominate  since  they  are 
important  near  the  source  while  the  long-range  contri¬ 
butions  of  the  discrete  modes  are  severely  attenuated. 

IV.  SUMMARY 

We  have  presented  a  model  of  surface  generated  noise 
in  the  ocean  in  which  the  random  noise  sources  are  rep¬ 
resented  by  correlated  monopoles  distributed  over  an 
infinite  plane  parallel  to,  and  located  on  arbitrary  depth 
below,  the  ocean  surface.  Expressions  have  been  de¬ 
rived  for  the  intensity  and  spatial  coherence  of  the  noise 
Reid  in  a  stratified  medium  based  on  a  wave-theoretical 
treatment.  Examples  have  been  given  which  demon¬ 
strate  that  environmental  factors,  such  as  the  sound- 
speed  profile  and  the  presence  of  the  bottom,  can  be 
important  in  determining  the  spatial  properties  of  the 
noise  field.  We  have  also  shown  that  for  an  isovelocity, 
semi-infinite  fluid  medium  our  results  are  identical  to 
those  of  previous  investigators. 
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APPENDIX  A:  SPATIAL  COHERENCE  IN  A 
HOMOGENEOUS  SEMI-INFINITE  SPACE 

In  this  appendix  we  show  analytically  that  the  theory 
developed  in  Sec.  I  reduces  to  earlier  work4  where  the 
ocean  was  modeled  as  a  semi-infinite  space.  For  this 
problem,  the  Green’s  function  is 

1  e*“  1  e***' 

G(r„r';z„z')  =  ^  -p-  .  <A1) 

where 

R=|  |r,-r'!* +  (*,-*')*]**,  (A2) 

and 


fl'  =  f  |r1-r'|*  +  U1+z'),),/* . 

and  k  -  ta/c,  where  o>  is  the  angular  frequency  of  the 
source  and  c  is  the  speed  of  sound  in  the  medium.  The 
Green’s  function  in  Eq.  (Al)  is  just  that  of  a  point 
source  and  its  image,  with  the  negative  sign  to  satisfy 
the  boundary  condition  at  the  surface, 

G(r1,r';0.i')  =  0.  (A3) 

To  obtain  g(tr,it,z')  we  note  that1* 

J_  e  ™ 

4s  R 

=  *jj~2  j  jfexpfilif-trj-r'l+qJzj-z'l)}-^  . 

<A4) 


where 

=  tOTk2>rf, 

~Hif-k2?a  for  k1  <if  . 

Thus  it  follows  from  Eq.  (8)  that 

=  i  /expttqJz.-zMl-expfaJri+r’m  ^ 

4b  \  tj,  / 

Similarly  we  have  that 

£*bj;z„z') 

4»  \  nl  /  ’ 

We  are  concerned  with  the  case  Zj.Zj  >z'.  Then,  from 
Eqs.  (A5)  and  (A6). 


1  exp[j(n,z,  -  n*z-)l  sin(n  z')sin(njz') 


.  (AT) 


Inserting  Eq.  (A7)  in  Eq.  (12)  we  obtain  the  expression 
far  the  cross-spectral  density  function 

Ci.tR.Zj.Zj) 

=  2s*cl  jf"  ^Cp)(f  q*jJo0j|R-p|) 

>-  e«Pf»(T,r.-n:.r«)lsin(n^,)sinfai;z,)\  ^ 

l»J,l  / 

If  k  is  real,  the  integral  inside  the  large  parentheses 
can  be  written 

f*  J0(ij  (R  -  p|  )exp(i(z,  —  Zj)(fe*  —  ij*)1*  1 


b  -  T] 

* /*  JofqlR-phexpf-Uj+ZjKi^-b1)’'^) 

* 

„  sinlr [z'(jf  -  b*)1'^  ]  ,  _ 

x - ^_ki - lndn-  tA9) 

We  now  must  make  some  assumption  about  S(  p). 

Cron  and  Sherman4  have  calculated  the  correlation  func- 
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N(p)=2mrnl(ltp)-aJm(kp) 


(A  10) 


tion  for  a  homogeneous  half-space  with  noise  sources 
having  cos'd  directionality  distributed  uniformly  over 
the  surface.  Liggett  and  Jacobson*  have  shown  that 
cos'd  directionality  is  equivalent  to  assuming  omnidi¬ 
rectional  sources  with  a  correlation  function  given  by 


for  m%l.  Inserting Eq.  (A  10)  in  Eq.  (A8),  the  angular 
part  of  the  p  integration  can  be  performed,  resulting  in 
the  expression 


2'mlfe-'^  F.(q)J0OjR)expf  iZ(k*  -  jfY'2] 


sittUW-r??*] 

k*-rf 


(  Fm(n)J0(j}R)expl-(z1+Zi)(i?-k2Y/x} 


sinli1  [  t  '{rf  — 


where  Z  =  z1-zs.  Fm(ij)  is  known,  and  is  given  by1 


/„(0 ,Z)=  f*  S**'1  co S(ZC)J£. 


Falrj)  =  f  Jm<£-p)J0(np)pi"'llP 
J0 

j  0  for  k<q, 

”  U2-,(*i-n2)]"V'[ri 


j  0  for  k<q, 

l  ^2-,(*i-7rt)]"V[^(M)^*,  for*>j,.  CA12) 
Thus,  using  Eq.  (A12),  Eq.  (All)  becomes 


Cu{R\el,zx)sq*2mk~*m  f  (**- if)"'lJ0(tiR) 
Jo 


x  expjiZ(t2  - 1]2  Y** ) 
„  sixtl z'(k3-rfY*) 


k*-rf 


For  m  =  1,  the  integral  is  easily  calculated.  The  result 


/t{0.Z)-kZ~1  sinGtiZ)  +  2T”1! coslfcZ)  —  1  ] .  (A21) 

Thus,  we  have,  for  m  =  l. 

C_.(0,Z)  =2(t2)'1  sin(feZ)  +  2<*2)*2[cos0t2)-  1  ] .  (A22) 

For  m  >1  ire  note  that 

IJ0.Z)  =  (-ir-1  -|^/1<0.2);  CA23) 


To  compare  with  Cron  and  Sherman’s  results,  we  let 
z'- 0  and  take  the  normalized  function  C^(R:z,.z~h 
Thus  (Re  denotes  the  real  part) 

CjR;z„  a*) 

‘4i4’ 

Then  we  must  evaluate  the  integral 

TjR,Z)-J  (^-lfr^JobiRicosIZ^-ifr^hdri. 


First  consider  the  case  when  z ,  =%.  Then  we  have 
fjR,0)=  (A  16) 

This  is  a  standard  integral  and  is  given  by” 

fm(R,0)=2m~lkmR~m(rrt—  1 )!«/_{«).  (A17) 

Therefore, 

CjR;z„z,)  =  2amlJm(kR)/(kRr,  CA18) 

which  is  just  the  correlation  function  of  the  surface,  Eq. 
(A  10). 

Next,  let  R=0.  Then  Eq.  (A  15)  reduces  to 

/.(0,Z)=  r  tf-Tfr-'cosm'-ifntb).  (A19) 

which,  after  changing  to  the  variable  £.  where  J 
=(feJ-ifJ)1'1,  becomes 
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so  for  m  >  1 


Ct,(0.2)  =  (-l)'-,2.fr-i 


-S' -2 

2* -I  f|t0,2)  . 


The  results  expressed  by  Eqs.  (A  18),  (A22),  and  (A24) 
are  in  agreement  with  those  of  Cron  and  Sherman. 

A  comment  is  in  order  about  the  significance  of  the 
second  integral  in  the  large  parentheses  of  Eq.  (All). 
Because  of  the  source  correlation  function  chosen  ( see 
Eq.  (A  10)].  the  second  integral  vanishes.  We  could 
have  taken  the  sources  to  be  completely  uncorrelated  by 
using  the  .V{p)  given  by  Eq.  (17);  then  each  source  would 
be  equivalent  to  an  independent  dipole.  The  first  term 
of  Eq.  (All)  would  then  give  Cron  and  Sherman’s  re¬ 
sults,  while  the  second  term,  which  can  easily  be  cal¬ 
culated,  would  be  negligible  except  near  the  surface. 
The  second  term  therefore  is  the  contribution  of  the  * 
nearfield  of  the  dipoles. 

APPENDIX  B:  A  SIMPLE  NORMAL-MODE  EXAMPLE 

As  an  illustration  of  the  normal-mode  representation 
of  the  noise  field  we  consider  an  isovelocity  waveguide 
of  depth  //  bounded  above  by  a  pressure- release  sur¬ 
face  and  below  by  a  rigid  bottom.  The  bewndary  condi¬ 
tions  for  this  problem  are 

C,(0)=0  (Bla) 


=0,  (Bib) 

I«*» 

where  the  functions  UJLz)  satisfy  Eq.  1 23).  The  solution 
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of  Eqs.  (23)  and  (Bl)  is 

UK(z)  =  {, 2/f/)l/asin(X^),  (B2) 

With 

k2n  =  k2-\2„,  (B3) 

where 

».-(aTi  i)'  ”-w .  ®4) 

We  introduce  attenuation  into  the  system  by  letting  the 
wavenumber  k  =  u/c  be  complex: 

k  =  K  +  te  (B5) 

with  e,  the  plane-wave  attenuation  coefficient,  taken  to 
be  a  small  positive  number.  Equation  (B3)  indicates 
that  the  modal  wavenumber  k„  must  also  be  complex. 
Thus,  we  let 

k*  =  K,  +  i<*n-  (B6> 

The  modal  attenuation  coefficient  or,  can  be  shown  by  a 
similar  method  to  that  of  the  Appendix  of  Ref.  15,  to  be 
given  by 

«„=—/  K[z)\Vn(z)\2dz,  (B7) 

Kn  */o 

for  the  general  case  of  a  depth-dependent  sound  velocity 
profile  c(z)  =  <<■/«.  In  the  case  considered  here,  c(z)  =  c 
a  constant,  so  Eq.  (B7)  reduces  to 

an  =  £k/k„.  (B8) 

Substituting  (B2)  and  (B8)  into  Eq.  (30)  we  obtain  the 
simple  result 

C„(R,  Zv  z2)=Zirq2p2sW)/ei<?H2  £  anJ(\„2')  sinU^) 

n 

*sin(A  „z2)J0{k„R).  (B9) 


We  note  that  the  attenuation  coefficient  c  appears  in 
Eq.  (B9)  just  as  a  scaling  factor  and  does  not  affect  the 
form  c*  the  cross-spectral  density  function.  This  is  a 
result  of  the  simple  example  chosen.  In  general,  the 
factor  (a,,*,)"1,  which  weighs  each  term  in  the  sum  will 
not  be  a  constant,  but  will  depend  on  n. 
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